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EFFECT OF ARTIFICIAL PITCH DAMPING ON THE LONGITUDINAL 
AND ROLLING STABILITY OF AIRCRAFT WITH 
NEGATIVE STATIC MARGINS 
By Martin T. Moul and Lawrence W. Brown .. 

SUMMARY 


A preliminary theoretical investigation has been made of the short- 
period longitudinal and steady-rolling (inertia coupling) stability of 
a hypersonic glider configuration for center-of -gravity locations rear- 
ward of the airplane neutral point. Such center-of-gravity positions 
for subsonic flight would improve performance by reducing supersonic 
and hypersonic static margins and trim drag. Results are presented of 
stability calculations and a simulator study for a velocity of 700 ft/sec 
and an altitude of 40,000 feet. 

With no augmentation, the airplane was rapidly divergent and was 
considered unsatisfactory in the simulator study. When a pitch damper 
was employed as a stability augmenter, the short-period mode became 
overdamped, and the airplane was easily controlled on the simulator. 

A steady-rolling analysis showed that the airplane can be made free of 
rolling divergence for all roll rates with an appropriate damper gain. 


INTRODUCTION 


Aircraft performance at supersonic and hypersonic speeds is some- 
times penalized by high trim drag. This unfavorable condition results 
from the rearward shift of the neutral point (stability increase) for 
supersonic speeds and the requirement of static stability at subsonic 
speeds . The possibility occurs of improving the performance of air- 
craft which are operated predominantly at high speeds by permitting 
some static instability at subsonic speeds. 
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As the center of gravity is moved rearward, the occurrence of prob- 
lems in both longitudinal stability and roll coupling can be anticipated 
at subsonic speeds. In the longitudinal mode, a pitch divergence is 
encountered when the center of gravity is located rearward of the maneu- 
ver point. In the rolling mode, it has been shown analytically (ref. 1) 
and demonstrated in flight (unintentionally in most instances) that roll 
coupling and a severe divergence may be encountered when the rolling 
velocity exceeds a certain critical value dependent on the pitching- 
moment-curve slope (static margin) or the yawing-moment-curve slope. 

Both of these stability problems must be investigated in any considera- 
tion of flight at small or negative static margins. 

A preliminary analysis has been made of the stability both in pitch 
and steady roll (inertia coupling) of a hypersonic glider configuration 
for a subsonic-speed condition in which the airplane center of gravity 
is rearward of the neutral point. Since most high-speed aircraft are 
equipped with pitch dampers, the effect of auxiliary pitch damping in 
stabilizing such a system is considered. 

Results of response calculations and a simulator study for a veloc- 
ity of 700 ft/sec and an altitude of L0,000 feet are included. 


SYMBOLS 


a n 

A 1 > B 1> C 1 

b 

B, C 


C L 


j m 


normal acceleration, g units 

coefficients of critical-rolling- velocity equation 
wing span, ft 

coefficients of longitudinal ciaracteristic equation 

wing mean aerodynamic chord, ft 

Lift 


lift coefficient, 


qS 


. . , . , . , Pitching moment 

pitching-moment coefficient, qQ~ 

. . , Yawing moment 

yawing-moment coefficient, 


lateral-force coefficient. 


Lateral force 
qS 



CCL 
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g 

h 

ix 

% 

iz 

K 

k y 

m 

p 0 

Per 

q. 

r 

S 

t 

t l/2 

V 

a 

d 


acceleration due to gravity, ft/sec^ 
altitude, ft 

moment of inertia of airplane about principal X-axis, 
slug-ft^ 

moment of inertia of airplane about principal Y-axis, 
slug-ft^ 

moment of inertia of airplane about principal Z-axis, 
slug-ft 2 

pitch-damper gain, deg/deg/sec 
nondimensional radius of gyration about Y-axis 

airplane mass, slugs 

steady rolling velocity, radians/sec 

critical rolling velocity, radians/sec 

dynamic pressure, lb/sq ft 
yawing velocity, radians/sec 
wing area, sq ft 
time, sec 

time to damp to one-half amplitude, sec 
airplane velocity, ft/sec 

angle of attack of airplane principal axis, radians or deg 

time rate of change of angle of attack, — , radians/sec 

dt 

initial angle of attack, deg 

angle of sideslip, radians 
elevator deflection, radians or deg 
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e,p 


A 

M- 


m 

pSc 


pilot-applied, elevator deflection, radians or deg 


of pitch, radian 

d 0 

radians/sec 


angle ui jjj. i etu.-Lciii 

pitching angular velocity, — , 

dt 

pitching angular acceleration, radians/sec 2 

dt 2 

characteristic root of quadratic equation 



SCl 

° L a 

ckx 

C i%, 

bet 

C m 


m q 



d 2V 

C n 

n p 

(5C n 
" c>p 

C n 

“r 

_ dC lx 

a rb 


d 2V 

Cv 

dCy 
' ^ 


atmospheric density , slugs/cu ft 


t = 


m 


pSV 


static margin 


C™ = C m + KC 
m le ra q c 


% 
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d&e 


L a - 1 SC L, 


'a 


= qScC 


% 


M - p 

<1 - 2V ^q 


N p = qSbC n , 


P 


Nr. = 


qSb^ 

2V 


Yp = qSC Y 


= M q + K^e 

»%e " 1 S5 S e 


ANALYSIS 


An airplane having a center-of-gravity position at or rearward of 
the neutral point and no stability augmentation may possess undesirable 
stability characteristics in both the pitch and roll modes. The effect 
of artificial pitch damping on the stability of these two modes will be 

considered . 


Longitudinal Stability 

As the short-period mode is the longitudinal mode of primary con- 
cern to pilots , this analysis has been restricted to a consideration of 
this mode. The airplane characteristics which influence pilots' opin- 
ions are period and damping of the oscillation and stick force per g. 
This investigation is limited to a consideration of the dynamic charac- 
teristics for variations in the aerodynamic parameters, pitching-moment- 
curve slope and damping in pitch. In the simulator study, satisfactory 
stick forces were mechanized so that pilots’ opinions were dependent 
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only on the airplane dynamic characteristics simulated. The dynamic 
characteristics are analyzed for zero and unstable (positive) values of 
the pitching-moment-curve slope with and without a pitch damper included 
as a stability augmenter. 


Equations of motion and system transfer function .- The longitudinal 
short-period characteristics are determined from the pitching-moment and 
lift equations of motion. These equations are written as 

I Y 0 = M^e + Ma,a + ^ e & e 
mV a - mV6 = -L a a 

when weight components, lift due to elevator deflection, pitching moment 
due to &, and variations in forward speed are neglected. 


a) 


To increase pitch damping, an autopilot introducing an elevator 
deflection proportional to pitching velocity (6 e = K0) is considered. 

Preliminary calculations and results of reference 2 showed that small 
gyro time lags have a negligible effect on airplane response. Conse- 
quently the effects of time lag are omitted in this analysis. The equa- 
tions of motion with artificial pitch damping are 


where 


I Y 0 = MqgG + Maa + M& e 6 e ,p 
mVa - mVG = -L a a 


(la) 


M % = + K^ e 


( 2 ) 


and 5^ ~ is a pilot-applied elevator deflection. The airplane trans- 
fer functions a/& e p and Q/S e p are 


a 


M8 e /3 


1 


e >P t. 2 , / Lq _ hit _ . LlS; 


^mV ly J \l\ ' mV ly 




M&e > . L a Mi»e 
Iv mV I-r 


e >P >2 , , L a ^3eV /Ht , ■'■’a \ 


A + : -S _ _ 

\mV l y 


i + 

I- mV Iv 


•' i 


(3: 



7 


Characteristic equation and stability .- The characteristic equation 
( denominator of the transfer function) is written generally as 

A 2 + BA + C = 0 (k) 

and has the roots 



For stability, the coefficients of a quadratic characteristic equation 
must be positive. Now, the coefficient B is always positive for nor- 
mal values of and C^, that is, positive Cj^ and negative 

values. The coefficient C is normally positive (for negative 

pitching-moment- curve slopes) but may become negative for positive values 
( ^m a * The condition of interest for this investigation is that of 

positive values, for which artificial pitch damping is employed 

to maintain positive values of coefficient C. This condition corre- 
sponds to a center- of- gravity location between the neutral point and 
maneuver point. The maneuver point is the center-of-gravity position 
for neutral maneuvering stability and is defined by C = 0. 

The nature of the characteristic modes is determined by the sign 
of" the discriminant (B/2) — C. Substituting for B and C gives 




L q, ^3e 
mV I Y 



(5) 


2 

must be positive or zero. Then for > 0, 

ive and the characteristic roots are real; and two 
aperiodic modes (which are stable for C > 0) rather than the usual 
short-period oscillation are obtained. These conditions > 0 and 

C > o) hold whenever the center of gravity is between the neutral and 
maneuver points and are the primary conditions of concern in this 
analysis . 


The term 


■ mV I y ) 
C is posit 
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Stability in Steady Roll 


Preliminary analyses of roll-coupling problems are frequently made 
by the use of the equations of motion for steady rolling and employ the 
concepts of critical roll rates, rolling divergence boundaries, and 
stability and airplane responses in steady rolls. (See refs. 1, 3, 
and 4.) In reference 1 it was shown that the rolling velocity, above 
which a divergence can be expected with zero values of damping, is given 


by the smaller value of 



For the conditions of 


interest in this paper, the longitudinal response consists of two aperi- 
odic modes by virtue of artificial pitch damping, and it is of interest 
to investigate the effect of rolling velocity on the stability of these 
configurations . 


The following constant coefficient of the characteristic equation 
for steady rolling ^including Mq^, N r , L a , Yp) is from reference k: 



M 3e N r o L a Yp fl x - I Y \ / l Z - Ixl p % Np 

I Y I z mV mV y I z ){ I Y r° “ Iy I Z 


Np L a N r Yp N r L a 

Iy 1 2 , mV Iy Iz, mV Iy IZj mV mV 


The coefficient E equated to zero will yield critical rolling veloc- 
ities as functions of aerodynamics and inertia characteristics. In 
general form, the equation is 


A lPcr^ + B lPcr 2 + Cp = 0 


from which 



Positive values of P cr ^ are indicative of regions 
bility, whereas negative or complex values of P cr ^ 

are indicative of configurations free of rolling divergence for all 
rolling velocities. 


of rolling insta- 
ller imaginary) 



9 


For stability-boundary considerations, there are two conditions of 
interest : 


(l) For negative Bp, the boundary between real and imaginary roll 
rates is defined by the discriminant equated to zero. Thus, 


£l 

Al 



( 8 ) 


(2) For positive Bp, the boundary between real and imaginary roll 
rates takes a different form; that is, p cr ^ ~ from which the result 
is obtained 



(9) 


With equations (8) and (9), stability boundaries in terms of selected 
aerodynamic derivatives as variables may be constructed. An application 
of these results will be given subsequently. 


RESULTS 


Calculations of longitudinal and rolling stability were made for 
a high-speed aircraft conf iguration for a velocity of 700 ft/sec and 
an altitude of 40,000 feet. Physical and aerodynamic characteristics 
assumed for this investigation are given in table I. Results are pre- 
sented for a range of static margins from 0.02 to -0.02 and pitch- 
damper gains from 0 to 2. 


Longitudinal Stability 

For positive damping conditions, there are three possible types of 
short-period modes: a damped oscillation, an overdamped response (two 

stable aperiodic modes), or a divergent and convergent mode, depending 
on the center-of-gravity position. In this investigation the center-of- 
gravity positions of interest are those falling up to 0.02c rearward of 
the neutral and unaugmented maneuver points, for which the characteris- 
tic motion with the addition of artificial pitch damping is an over- 
damped response. 

Constant t curves.- Time to damp to one-half amplitude is the 
characteristic used to define overdamped systems . An equation relating 
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determined by substituting A = - 2 

x '\/2 

into the characteristic equation (eq. (4)}. The result is 


t l/2 to C m qe and Sx iS 


2 > xt l /2 


+ 1 -? 86Tt l/2 K Y 2c L a ~ 1-92% 2 t 2 


\ o .693 t - | C La t l/2 l 0.693Tt l/2 - | c La t l/2 


Lines of constant t ^/2 are P resen ted in figure 1 with as 

the ordinate amd as the abscissa. There are two characteristic 

% 

damping times associated with each point as indicated by the 

intersecting lines. 


Within the parabola in the right-hand quadrant are contained the 
conditions for an oscillatory short-period response. The divergence 
boundary, or maneuver-point location, appears in the left-hand quadrant. 
Between the oscillatory and divergence boundaries is the region for 
stable overdamped responses. Stability deteriorates as static insta- 
bility ( positive values j increases or decreases, as evi- 


v 

denced by the larger t 


1/2 


values . 


Effect of damper gain on twg •" figure 2 is presented the varia- 

L 

tion of tp /2 w ith damper gain K for values of — 21 0 f 0, 0.01, 


dCn 


and 0.02. The upper set of curves for — ^ corresponds to the dominant 

dC-L 

mode in the longitudinal response, and it is desirable to keep ^\j2. 
small. As gain approaches infinity, the upper set of curves becomes 


asymptotic to the value 0.693 


mV 


a 


For this configuration, there is 


little improvement in damping for increases in gain above 3- 

The effect of artificial damping on the controllability of this 
particular configuration is indicated in the figure by the curves of 

cL Cw> 

constant a/o e . Both pitch damping and — — contribute to the airplane 

c:C L 

stiffness (coefficient C) so that the magnitude of a/6 e decreases 


with an increase in K or a decrease in 


ac m 

dCr ‘ 


The usual range of a/& 6 


for airplanes is about -1 to -3. Thus, figure 2 is useful for selecting 
pitch-damper gains in accordance with any selected stability and control 
requirements . 
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Response to unit-step elevator deflection .- The dynamics of several 
configurations are illustrated in figure 3 by means of angle-of-attack 
responses to unit-step elevator deflections for several static margins 
and pitch-damper gains. With no artificial damping, the maneuver point 
and neutral point are nearly coincident. Hence, the response for K = 0 


is divergent for = 0.02, about neutrally stable for = 0, and 


dC L 

stable and aerodynamically damped for 


ac m 

dCr 


dC L 

= -0.02. With artificial 


pitch damping added to the divergent configurations, the responses 
become stable and overdamped, as shown in figure 1, and would be expected 
to be tolerable to pilots . 


Simulator results .- In order to determine a pilot* s opinion of 
negative static-margin configurations both with and without artificial 
damping a simple closed- loop simulator consisting of an analog computer, 
chair, center control stick, and displays of normal acceleration and 
angle of attack was utilized. The stick was mechanized to provide a 

force per g of 4^ pounds and a deflection per g of 1 inch for stable 

^ dC 

conf igurations . For the cases of no artificial damping and — — - o 

dC L 

and 0.02, the stick force per g was about 0 and -3 pounds (push), 
respectively. The pilot’s task was to pull up from 1 g to 3g and trim 
the airplane. This task and the simulation of the airplane by the two- 
degree-of-freedom equations of motion are consistent and adequate for 
an analysis of short-period dynamics. 


The simulator results, as time histories of normal acceleration 
and pilot’s input to the elevator, are presented in figure 4 for two 


static margins, 


dC- 


m 


dC L 


= 0 and 0.02. In figure 4(a), the results for the 


dC n 


basic airplane with no artificial damping are shown. With — — = 0 

dC L 

dCm 

the airplane Is about neutrally stable, and with — ^ = 0.02, the air- 

dC L 

plane is rapidly divergent. For both configurations the pilot was able 
to pull up to 3g and steady out. 


With a static margin of 0, no pitching moment is produced by angle 
of attack. Consequently, the stick was first pulled back to pitch the 
airplane and then almost neutralized when the desired normal accelera- 
tion was obtained. The pilot was able to perform this task without 
difficulty although the control motions and forces are different from 
those required for maneuvering an airplane with positive static margins 
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With negative static margin the task was difficult and the record 
was taken only after several practice runs . After initiating the motion 
with the control, the pilot had to rapidly reverse the elevator deflec- 
tion in order to halt the divergence and trim the nose-up pitching 
moment. Although he could control the motion, the pilot considered this 
configuration unflyable for operational conditions because the control 
task required about 90 percent of his time. 

With artificial pitch damping added, the unit-step responses for 

dCm 

— 2 = 0 and 0.02 are shown in figure 3 and the simulator results are 
dC L 

shown in figure 4(b) . The pilot found that these systems were easy to 
control. Since these configurations were overdamped, the pilot could 
pull up rapidly by initially overcontrolling and then easing off on the 
control as the desired normal-acceleration level was approached. In 

figure 4(b), the response for = 0.02 is slower because the pilot 

dCp 

did not overcontrol as much as for the — S = 0 condition. In this 

dC L 

application of a pitch damper, customary washout circuits must be omitted 
in order that the damper remain effective in the steady state. As a 
result of this simulator study, it appears that hypersonic glider con- 
figurations may be provided with satisfactory short-period characteris- 
tics by a pitch damper for center- of -grav: ty positions behind the neutral 
point . 


Roll Coupling 


Stability in steady roll .- Conditions were given previously for 
avoiding rolling divergence over the entiie range of rolling velocities. 
(See eqs . (8) and (9)*) With the given relations, boundary curves have 
been determined as a function of pitch-damper gain and pitching-moment- 
curve slope, the two significant parameters in this investigation, and 
are presented in figure 5* For values of K above the boundary curve 
for rolling stability, the particular configuration will be free of 
rolling divergence. The corresponding boindary for neutral longitudinal 
stability (lower curve) is also shown. A1 small values of C^, r °Fling 

stability requires a somewhat higher pitch-damper gain than does longi- 
tudinal stability. 


The region between the curves corresponds to conditions of C m 

lu a 

and K for which rolling divergence would be expected for some finite 
range of roll rates. Rolling instabilities of this kind are character- 
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istic of many supersonic fighter aircraft and analyses are usually made 
of specific conditions by means of rolling divergence boundary plots. 

For examples, see references 1 , 3 , and 4. 

To further illustrate the characteristics of the region between 
the two boundary curves of figure 5 * critical rolling velocities are 
presented in figure 6 as a function of a damper gain for = 0 , 

0.015, and 0.030. The curves originate at the p cr = 0 axis at the 

minimum K value for longitudinal stability given in figure 5 • Inter- 
sections of a given curve with a particular value of K define 

the roll-rate range for which a specific configuration would be diver- 
gent . The upper portions of the curves extend to K values less than 
the minimum for longitudinal stability and indicate the condition of 
spin stabilization of longitudinally divergent configurations. 

The characteristics roots for roll rates of 0.5, 1.0, and 
2.0 radians/sec are given in table II for four conditions of C m 
and K. 

The configurations having C m = 0.04, K = 1.88 and C =0, 

K = 0.48 would appear in the completely stable zone of figure 5, and 
the characteristic roots given in the table are all stable. The other 
two configurations having = 0.04, K = 0.83 and = 0, K = 0.15 

would appear in the area between the boundaries in figure 5 and each 
would have a divergent mode for a range of roll rates . In table II a 
divergent mode is indicated for each of these configurations at a 
particular roll rate. 

Transient responses .- Two examples of angle-of-attack responses 
during steady rolling are presented in figure 7 for a roll rate of 
0.5 radian/sec and an initial angle of attack of 5 °. The configurations 
selected were = 0.04, K = 0 . 83 and = 0.04, K = 1.88 for 

which the characteristic roots are given in table II. For K = 1.88 
the response is stable with small amplitudes. For K = O. 83 , the 
response is slowly divergent but may not be objectionable to pilots. 


Phugoid Characteristics 

The short-period oscillation is usually the longitudinal mode of 
primary concern to pilots and the purpose of this investigation was to 
investigate a method of improving its characteristics. Although pitch 
damping stabilized the short-period mode, examination of the quartic 
longitudinal stability equation, obtained when perturbations in forward 
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speed are permitted, indicates a divergent condition for unstable 
pitching-moment curves and a negligible effect of pitch damping on 
this instability. This divergence is associated with the phugoid char- 
acteristics and flight tests of such configurations are required to 
determine the significance of this divergence. If the phugoid charac- 
teristics prove to be undesirable, an additional feedback of a quantity 
significant in the phugoid mode, such as forward speed, is required. 

An alternate control system, employing normal acceleration or angle- of- 
attack feedback, may be employed to improve simultaneously both short- 
period and phugoid characteristics. 

SUMMARY OF RESULTS 


A brief analysis, by means of stability calculations and a simulator 
study, has been made of the stability and controllability of a hypersonic 
glider configuration for center-of-gravity locations rearward of the 
neutral point. When a pitch damper was employed as a stability augmenter, 
the following results were obtained for a velocity of 700 ft/sec and an 
altitude of 40,000 feet: 

1. Without augmentation, the airplane was rapidly divergent for 
negative static margins and was considered unsatisfactory in the simu- 
lator study. 

2. With artificial pitch damping, the short-period response is 
stable and overdamped and the configuration was easily controlled on 
the simulator. An instability of the phugoid mode exists, however, for 
negative static margins and investigation by simulator or flight studies 
is required to determine the effect of this mode on the controllability 
of the airplane . 

3. A roll-coupling analysis showed that the airplane can be made 
stable for all roll rates with an appropriate damper gain. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., February 5, 1959- 
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TABLE I 

AIRCRAFT PHYSICAL AND AERODYNAMIC 
CHARACTERISTICS 


in, slugs 585 

I X , slug- ft 2 11,000 

l y, SlUg- ft 2 126,000 

l z , slug- ft 2 136,000 

b, ft 35 

c, ft 25 

S, sq ft 750 

h, ft 40,000 

V, ft/sec 700 

Cp^, per radian 2 

Cm^, per radian -0.6 

C^ , per radian -0.08 

Cy^ , per radian -0 . 286 

C n p, per radian 0.057 

C nr , per radian -0.31 



TABLE II 


CHARACTERISTIC ROOTS IN STEADY ROLLING 


0.04 1.88 


0.48 


Po> 

radians/sec 


Roots 

-0.224, 

-3.645, 

-0.186 ± 

1.3581 

-0.236, 

-3-489, 

-0.260 ± 

1.6401 

-0.931, 

-2.591, 

-0.361 ± 

2.438i 

+0.013, 

-2.095, 

-0.184 ± 

1.3791 

-0.072, 

-1.874, 

-0.253 ± 

1.6951 

-0.902 d 

0.836i, 

-O.323 1 

2.5131 

-0.245, 

-I.151, 

-O.232 ± 

1.4371 

-0.103, 

-1.124, 

-O.316 ± 

1.8541 

-0.555 i 

: 0.9131, 

-0.375 l 

: 2.7991 

-0.054, 

-0.840, 

-0.201 ± 

1.4641 

+0.058, 

-0.858, 

-0.249 ± 

1.8871 

-O.366 i 

: 0.9361, 

-0.282 ± 

2 . 826i 
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3.- Angle -of-attack response to unit-step elevator deflection. 





tjSec. 

(a) No artificial damping. 

Figure 4.- Simulator results for a prescribed control task with rear 

ward center-of-gravity positions. 









